ADAOyiSO? 


DEPARTMENT  OF  DEFENCE 

AUSTRALIAN  DEFENCE  SCIENTIFIC  SERVICE 
MATERIALS  RESEARCH  LABORATORIES  ^ 

MARIBYRNONG  VICTORIA 

TECHNICAL  NOTE 


MRL-TN-371 


/ 


FORM-FUNCTIONS  FOR  MULTI-COMPONENT  PROPELLANT  CHARGES 
INCLUDING  INHIBITED  GRAINS  AND  SLIVER  BURN 


v/ 


J.  Stals 


APPRO  VI.  D 

n 

! 

pi  True  RPivur 

/-\L 

MELBOURNE 


d r 

r3!  r. 


^ NOV  8 19TB 


,HT 


vCO!  iiiw.iwedllli  of  Aui'ra'ia 

SEPTEMBER,  1975 


S1C1ED 


f 

1 H7  U\H’:o  J-TAT-. »:*•! 

i 

J 

T r.CHM’CAL  P.:  'MAT... M ' .. 

! 

;S  Aur®;  ‘,-0  10 

REPROL  UCi.  AN:'  5 ;.U.  Th1;;  i, 

'.u>V  J 

1 1 


1.  Security  Grading/Release 
Limitation  : 

(a)  Document  content  : 

Unclassified 

(b)  This  Page  : 

Unclassified 


2.  Document  Type /Number  : 

MRL-TN- 371 

3.  Document  Date  : 

SEPTEMBER,  1975 


»• 

i 


4.  Title  and  Sub-Title  : 

FORM-FUNCTIONS  FOR  MULTI -COMPONENT  PROPELLANT  CHARGES 
INCLUDING  INHIBITED  GRAINS  AND  SLIVER  BURN 


5.  Personal  Author 


J.  Stals 


6.  Corporate  Author 


Materials  Research  Laboratories, 
Maribyrnong,  Victoria 


7.  Abstract  (See  reverse) 


8.  Descriptors  : 


9.  Cosati  Classifications  : 


Propellant  Grains 
Grain  Shapes 
Area 

Inhibition 
Burning  Rate 
Interior  Ballistics 


1901 


11.  Task  Reference  (RD/P8 ) : 
TASK  71/11 


j 


10.  Library  distribution  : 

MRL  (Maribyrnong) 

MRL  (S. A.  Branch) 

MRL  (N.S.W.  Branch) 

WRE 

Russell  Library 
Campbell  Park  Library 


12.  Sponsoring  Agency 
Reference  : 


13.  Cost  Code  : 
21804 


f 


14.  Imprint 


MARIBYRNONG  - MRL  - 1975 


I 


ABSTRACT 


V 

Cubic  form- functions  in  burning  surface  regression  distance 
are  derived  and  compiled  for  granulations  that  are  used  in 
current  U.S.A. , U.K.  and  Australian  ammunition.  Functions  for 
specific  inhibition  of  grain  surface,  sliver  burn  for  multi- 
perforated  grains,  and  unequal  webs  in  seven  perforated  grains, 
are  included.  The  functions  are  of  general  use  in  propellant 
grain  design,  closed-chamber  and  interior  ballistic  computations. 
Their  resolution  into  igniter  and  propellant  segment  contribu- 
tions as  a function  of  grain  surface  regression,  facilitates 
studies  on  multi-granular,  multi-propellant  charges,  mortar 
charges,  deterred  grains,  programmed  or  non-isochronous  ignition, 
and  the  design  of  specific,  or  maximum  progressivity  granulations. 
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NOTATION 


) Geometric  grain  constants 


Differential  of  z with  respect  to  r 

Mass  of  propellant  segment  i 

Initial  volume  available  to  propellant  gases 

Outer  grain  diameter  or  width 

Perforation  diameter  or  smaller  diameter  of  grain 
Packing  efficiency  of  propellant  charge 
Fraction  of  grain  web  unburnt 
Total  number  of  propellant  segments 


Form- factors 


Slit  length  in  slotted-hole  flake  grains 
Length  of  grain 

Number  of  grains  in  propellant  segment  i 

Number  of  grain  perforations 

Differential  of  S with  respect  to  r 

Grain  surface  regression  distance 

Ratio  of  grain  surface  area  to  its  initial  value, 
i.e.  grain  progressivity 

Slit  width  in  slotted  hole  flake  grains 

Surface  area  that  is  burning  at  time  t 

Initial  surface  area  of  the  grain 


Time  interval  after  primer  ignition 
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Initial  volume  of  propellant  grain 
Volume  of  propellant  grain  at  time  t 
Free  volume  left  in  charge  container 
Grain  web 

Fraction  of  charge  burnt  at  time  t 
Weight  fraction  of  segment  i in  total  charge 
Grain  density 
Form-fac tor 

Subtended  angles  as  shown  in  grain  sliver 
burn  diagrams 

Constant 

Spheroid  eccentricity  factor 


Subscripts 


e,E 
El, EE  ) 


) Denote  grain  end  surfaces  in  slivers 


I 

1,L 

LI, LL  ) 
max 


Segment  i of  total  propulsion  charge 
Primer  or  igniter 


m 


m 


m 


kg  m 
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rad 


) 

) Denote  grain  lateral  surfaces  in  slivers 

Maximum  value  of  subscripted  variable 
i,o  Inner  and  outer  webs 

s Segment  variable 

w Variable  value  when  r = web/2 

1 for  r Variable  value  when  small  sliver  is  all  burnt 

2 for  r Variable  value  when  large  sliver  is  all  burnt 

3 for  r Variable  value  when  residual  slivers  are  all  burnt,  N = 19 
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FORM-FUNCTIONS  FOR  MULTI -COMPONENT  PROPELLANT  CHARGES 


INCLUDING  INHIBITED  GRAINS  AND  SLIVER  BURN 


A.  INTRODUCTION 


The  optimisation  of  the  gasification  rate  for  a propulsion  charge  and 
the  design  of  granulations  to  achieve  it,  and  the  desired  weapon  performance, 
requires  in  part,  a quantitative  formulation  of  the  total  propulsion  charge 
form- function.  Namely,  the  fraction  of  charge  burnt  (z)  at  any  time  (t) , 
as  a function  of  grain  geometry  and  surface  regression  distance. 

In  the  past,  z has  been  expressed  in  terms  of  : 

1 2 7-11  17 

(i)  The  fraction  of  web  unburnt  (f)  ' ' ' 


z = (1-f)  (l+0f)  (1) 

or  z = K0  + Kjd-f)  + K2(l-f)2  + K3(l-f)3  (2) 

where  the  form-factors  0,  KQ  to  were  approximated  from  grain  geometry  or 
the  empirical  fitting  of  experimental  pressure-time  records.  For  closed- 
chamber  data,  the  geometric  values  of  0 are  generally  0.2-0. 3 too  low  as  a 
result  of  grain  heterogeneities  and  the  violation  of  Piobert’s  law  (see 
Section  B) . In  a gun,  however,  erosive  burning  results  in  greater  than 
predicted  grain  surface  regression  rates  and  gives  a compensating  approxima- 
tion for  the  variation  of  z with  time. 

In  addition,  the  variable  f loses  physical  reality  for  multiperforated , 
and  other,  grain  shapes  that  result  in  slivers  upon  web  burnthrough.  The 
form-factors  0 and  KQ  to  are  adjusted  empirically  to  compensate  for  the 
latter  sliver  bum  phenomena,  or  a compensating  ballistic  web  size  is 
invoked. 

(ii)  The  surface  area  (S)  of  the  grain1 

z = ( 1- (S/S  )2)  (1  + 0) 2/40  (3) 

o 

where  S is  the  initial  surface  area  of  the  grain  and  0 is  a form-factor. 
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The  treatment  of  geometric  and  chemical  heterogeneities  of  the  grain, 
its  specific  inhibited  surfaces,  in  terms  of  surface  area,  is  inconvenient. 

(iii)  The  regression  distance  (r)  of  the  burning  grain  surface. ^ ^,12  16 
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where  the  grain  geometric  coefficients  a.  are  explicitly  expressed  in  terms 
of  grain  dimensions. 

Equality  (4)  has  physical  reality  for  all  grain  shapes  and  it  can  be 
readily  used  to  treat  radial  chemical  gradients  in  grains,  such  as  those  due 
to  deterrents  and  inhibitors.  The  regression  distance  may  be  measured 
experimentally  upon  extinguishment  of  the  grain  flame,  and  its  time  deriva- 
tive is  the  useful  grain  surface  regression  rate  (dr/dt) . In  addition, 
sliver  burn  may  be  explicitly  formulated  in  terms  of  r.  Further,  the  use  of 
the  explicit  geometric  constants  a.  facilitates  parametric  studies  of  the 
effect  of  grain  geometry  on  propulsion  charge  performance  in  the  weapon. 

Since  a cubic  polynomial  in  r is  sufficient  to  define  most  grain 
geometries,  we  shall  use  Eq  (5)  as  the  basic  form-function  in  this  work  and 
relate  it  to  other  grain  properties  in  Section  B. 
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3ir 


d2r 


a3r 


(5) 


The  geometric  constants  of  the  grain,  viz.  a^,  a^  and  a^,  are  derived 
and  listed  in  the  Appendix  for  the  following  grain  shapes.  Selective 
surface  inhibition  is  allowed  for  in  (b)  and  (d)-type  grain  shapes. 


(a) 

Solid  right  prismatic 

(i) 

rectangular  and  ribbon 

(ii) 

cubic 

(iii) 

square  flake 

(b) 

Mortar  propellants 

(i) 

disc  and  cord 

(ii) 

slotted-'nole  flake 

(c) 

Ellipsoids 

(i) 

sphere 

(ii) 

oblate  spheroid 

(iii) 

prolate  spheroid 

(iv) 

symmetric  rolled-ball 

(d) 

Cylindrical 

(i) 

cord , of  b (i) 

(ii) 

tubular 

(iii) 

multiperforated  (N)  for  : 

2 


I N = 0,1,2 


",  prior  to  sliver  formation  and  web  burnthrough. 


II  N = 7 , 1 9 


during  sliver  burn, 


II  N = 7 


unequal  webs,  prior  to  and  during  sliver  burn. 


(e)  Slotted  tube 

The  above  shapes  cover  the  granulations  that  are  used  in  Australian 
ammunition.  They  also  incorporate  most  USA  and  UK  gun  and  mortar  pro- 
pellant configurations.  Specific  rocket  motor  grain  configurations  have 
been  treated  in  Refs.  19  to  25.  A useful  handbook  of  equations  for  mass 
and  area  property  calculations  for  various  geometric  shapes  is  listed  in 
Ref.  27. 

Some  of  the  form  functions  for  the  more  common  granulations  have  been 
published  previously ; ^ they  are  included  for  completeness.  The 
compilation  of  functions  given  in  the  appendix  forms  a useful  reference  for 
R & D work  on  propulsion  charge  design,  interior  ballistic  and  closed 
chamber  computations,  and  eliminates  the  need  for  time  consuming  derivation 
thereof. 


FORM-FUNCTION  AND  CHARGE  PARAMETERS 


Historically,  it  has  been  assumed  that  a solid  propellant  grain  burns 

under  the  following  conditions ; 

(a)  Each  grain  is  chemically  and  physically  homogeneous  in  all  directions 
and  there  is  no  cracking  or  mechanical  degradation  of  the  grain  during 
the  ignition  and  combustion  processes. 

(b)  The  segment  of  the  total  charge  in  question  consists  of  chemically  and 
geometrically  identical  grains  and  these  ionite  simultaneously  and 
uniformly  within  the  entire  segment  of  the  charge. 

(c)  Each  grain  is  ignited  instantaneously  over  all  its  uninhibited  external 
surface  (P.  Vieille,6  1893),  with  the  speed  of  igniter  flame  spread 
being  much  greater  than  the  initial  propellant  burning  rate. 

(d)  The  uninhibited  grain  surface  burns  in  parallel  layers  and  perpendi- 
cular to  the  surface  with  the  same  surface  regression  rate  over  the 
entire  burning  surface  (G.  Piobert,^  1839). 


Each  of  the  above  assumptions  are  violated  to  varying  degrees.  For 
example,  (a)  and  (b)  hold  under  ideal  propellant  production  conditions. 

These  are  rarely  achieved,  and  difficult  to  reproduce.  Vieille's  assumption 
and  (b)  hold  accurately  for  porous  propellant  beds  such  as  artillery  charges; 
they  break  down  however  for  compact,  high  loading  density  small  arms  charges 
in  which  non-isochronous  ignition  may  occur.  Piobert's  law  is  only  semi- 
quantitative,  particularly  within  perforations  of  multiperforated  grains, 
and  is  strongly  dependent  on  (a)  to  (c). 
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The  above  violations  cannot  be  easily  quantified,  nor  would  they  be 
statistically  reproducible  between  lots  of  the  same  propellant.  Assumptions 
(a)  to  (d)  are  thus  retained  and  used  to  provide  geometric  form  functions 
that  are  useful,  if  semi-quantitative,  in  propellant  grain  design,  the 
control  of  propellant  mass  burning  rate,  and  subsequent  computations  on 
interior  ballistics  and  closed  chamber  performance. 


The  fraction  of  charge  that  is  burnt  (z)  at  any  time  (t),  is  the  sum 


of  the  primer  (z^)  and  propellant  component  (z^)  contributions 


z = 5j.z 


I + Yj 


i = 1,2, 


(6) 


i=l 


where  the  total  charge  is  made  up  of  j segments  that  may  vary  in  chemical 
formulation  and/or  grain  geometry  and  dimensions.  The  weight  proportion  of 


igniter,  I,  and  propellant  component,  i,  are  denoted  by  6 and  6. 
respectively. 


If  V.  and  V . denote  the  respective  volumes  of  a grain  at  time  t and 

1 u 01 

mitially,  we  have 


z.  = 1.0  - V./V  . 
1 1 oi 


(7) 


and  for  segment  i of  the  propellant  charge,  Eq  (5)  becomes 


2 3 

z,  = an .r  + a„.r  + a..r 
l li  2i  3i 


(8) 


Upon  differentiating  Eq  (8)  with  respect  to  r we  obtain 


dzi/dr  = ali  + 2a2ir  + 3a3ir  E Bi  (S3y) 


(9) 


If  S . is  the  initial  surface  area  of  the  grain  and  S.  is  the  burning 
surface  area  at  time  t,  then 


a. . = S ./V  . 
li  oi  oi 


(10) 


and 


B.  = S./V  . 
1 l oi 


(11) 


The  ratio  (R^)  of  the  grain  surface  area  at  time  t to  its  initial 
value  is  thus 


R.  = S./S  . = B. /a, . 
1 1 01  1 li 


hence  a maximum  value  of  R.  exists  if  a,.  < o,  o r = - a„./3a  . r 

l 3i  2i  3i  max 


and  a,.  > o,  thus; 

2i 


R.  = 1.0  - a..  /3a,. a,, 
l max  2i  3i  li 


Let  us  define  P as  the  incremental  variation  of  S as  a function  of  r, 
that  is 


P . = dS./dr  = V . (2a, . + 6a,.r) 
ri  l oi  2i  3i 


P . = S . (2a  . + 6a  . r ) /a  . 
n oi  2i  3i  li 


where  for  a given  regression  distance  r the  grain  is  : 

(a)  Progressive  if  Pr  > 0.0  . S increases  with  increase  in  r 

to  give  an  increased  propellant  gasification  rate  and  burning 
surface  area. 

(b)  Neutral  for  P = 0.0  . A constant  burning  surface  and 
propellant  gasification  rate  is  maintained. 

(c)  Degressive  if  P < 0.0  . The  burning  surface  and  gas  evolution 
rate  decrease  with  increase  in  r and  t. 

If  p.,  C.  and  n,  denote  grain  density,  mass  and  number  of  grains 
respectiveiy  in  segment  i of  the  charge  then; 


n.  = C./p.  V . 
1 1 1 oi 


and  the  segment  charge  variables  (sub  s)  for  segment  i become 


S . = n.  S. 

SI  11 


P . = n.  P . 
sn  i n 


Let  Cap,  EP  and  denote  the  volume  available  (e.g.  chamber, 
cartridge,  etc.)  to  the  propellant  charge  prior  to  shot-start,  the  packing 
efficiency  and  free  air-space  respectively,  then 


V = Cap  - V 

Cl 


EP  = V /Cap 

d 


The  variables  a , a 2>  hence  B,  S,  R,  P^  and  z,  are  derived  in  the 

Appendix  for  the  grain  shapes  listed  in  Section  A. 


C.  CONCLUSIONS 


The  best  form  of  granulation  ballistically , is  the  one  that  will  give 
the  highest  muzzle  velocity  for  a given  charge  weight  and  chemical  composi- 
tion, within  the  tube  pressure  limits  for  which  the  weapon  was  designed. 

The  most  progressive  grain  shape  is  therefore  desirable.  The  listed  form- 
functions  permit  the  determination  of  the  granulation,  or  combinations 
thereof,  that  will  produce  a charge  of  a given,  or  maximum,  progressivity . 

In  addition,  the  compiled  form- functions  and  charge  parameters  facili- 
tate the  study  of  : 

(a)  Multi-granular  and  multi-propellant  charges,  e.g.  zoned 
artillery  ammunition,  blends  of  small  arms  powders. 

(b)  The  importance  of  the  contribution  of  the  primer  to  the 
propulsion  charge  energy,  e.g.  mortar  charges,  propellant 
actuated  devices. 

(c)  The  effect  of  inhibiting  specific  grain  surfaces,  e.g. 
rocket  motors,  propellant  actuated  devices,  very  progressive 
grains. 

(d)  The  presence  of  radial  chemical  gradients  within  the  grain, 
e.g.  deterred  small  arms  grains;  nitroglycerine  migration; 
solvent,  plasticizer  and  moisture  gradients;  sandwich  or 
layered  grains. 

(e)  Slivei  burn,  e.g.  multiperforated  grains. 

(f)  Unequal  webs  in  multiperforated  grains. 

(g)  Non- isochronous  charge  ignition,  e.g.  small  arms  charges, 
composite  and  moulded  charges,  programmed  ignition,  caseless 
ammunition. 


(h)  The  modular  use  of  the  charge  mass  burning  rate  law  in  charge 
design,  closed  chamber  and  interior  ballistic  simulations. 
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The  effects  of  non-uniform  erosive  burning  of  the  grain  surface  (e.g. 
within  perforations),  and  the  statistical  variations  in  grain  shape, 
dimensions,  perforation  asymmetry  and  chemical  composition,  have  been 
f omitted  from  the  compiled  form- functions  and  charge  parameters. 
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F.  APPENDICES 


Hence 

where 


2 3 

1-V/V  = a,r  + a_r  + a_,r 

o 1 2 3 

2 (LD+Lw+Dw) /LDw 
-4 (L+D+w) /LDw 


S =2  (LD  + Lw  + Dw) 
o 

V = LDw 

o 

S = 2C (L-u)  (D-u)  + (L-u)  (w-u)  + (D-u)  (w-u) 3 

V = (L-u) (D-u) (w-u) 

u = 2r 


(i)  Uninhibited 


s 

o 

= 

7i  D(O.SD+w) 

V 

o 

= 

2 

0.25  it  D w 

Hence 

S 

— 

7i  <D-u)  f 0.  5 (D-u)  + (w-u)  3 

V 

= 

0.  25  77  (D-u)2  (w-u) 

and 

z 

= 

1-V/V,  = + a^r2  + a3r3 

for 

al 

= 

2 (D+2w)/Dw 

a2 

= 

-4 (2D+w)/D2w 

a3 

= 

8/D2w 

(ii) 

Inhibited  Wall  v is  unchanaed 

S 

o 

= 

O 

0.5  Ti  D2 

Hence 

S 

= 

2 

0.5  7T  D - S i.e.  constant  burning 

surface  area 

V 

= 

0.  25  77  D2  (w-u) 

z 

= 

1-V/V  = a.r  + a.r2  + a r3 

ol2  3 

when 

32 

= 

a^  = 0.0 

31 

= 

2/w 

(iii)  Inhibited  Ends 


Vq  is  unchanged 


s 

o 

- 

TT  DW 

Hence 

S 

= 

r w(D-u) 

V 

= 

0.  25  tt  w(D-u) 

and 

z 

= 

i-V/Vo  = aLr 

when 

a3 

= 

0.0 

al 

= 

4/D 

a2 

= 

-4/D2 

(b) 

Slotted-Hole 

Flake 

w 


d + 21 


Assume  1 > w,  i.e.  no  slivers  are  formed. 


S = 2 (LD+Lw+Dw+lw-sw-ls)  + it  d(w-0.5d) 
o 

V = w (LD-  ls-O.  25  tt  d2) 

0 

Hence  S = 2C  L-u)  (D-u) - (1-u)  (s+u)  + (w-u)  (L+D+l-s-4  ) 1 + it  (d+u) 
■ [ (w-u)  -0.5  (d+u) 1 

V = (w-u)  C (L-u)  (D-u)  - (l-u)  (s+u)  - 0.25  it  (d+u)2l 

and  z = 1-V/V  = a.r  + a_r2  + a r3 

o 1 2 J 

for  a.  = 2w  (L+D+l-s+0. 5 tt  d)/V  + 2/w 

1 o 

a2  = -4  (2w+L+D+l-s-0. 25  tt  w+0.  5 tt  d)/VQ 


a3  = 2(8  - tt)Aq 
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Elevation 


Side  view 


2 2 *3 

C = (1-w  /D  ) 

S = 0.5  it  w(w  + D arcsin(t)/t) 

V = TT  Dw2/6 
o 

eu  = Cl- (w-u)2/ (D-u) ^ 

S =0.5  it  (w-u)  L (w-u)  + (D-u)  arcsin  (eu)/eul 

V = it  (D-u)  (w-u)  2/6 

2 3 

z = 1-V/V  = a.r  + a_r  + a_,r 

ol2  3 

= (4wD+2w2) /Dw2 
a^  = -(4D+8w)/Dw2 
a^  = 8/Dw2 

(d)  Rolled-Ball  : symmetric  sphere  of  radius  D/2  cut  by  two 

parallel  planes  at  a distance  w apart. 


2 

S = tt  Dw  + 0.5  i d 
o 

V = 0.25  TT  Wd2  + TT  w3/6 


Elevation 


Hence 


S = ti(D-u)(w-u)  +0.5  nLd-u  (w-D+d)/w3 


and 

for 


V = 0.25  it  (w-u)  Cd-U  (w-D+d)/w32  + 7!  (w-u)  ^/6 

2 3 

z = 1-V/V  = a,r  + a r + a r 

o 1 2 3 

a = [w2-d (D-w-1. 5d) 3/Cw (0.25d2+w2/6) 3 

a2  = C2D-3w-4d- (D-3d) (D-d)/wl/£w(0.25d2+w2/6) 1 

a3  = [5/3-2(D-d)/w+(D2-d2)/w23Aw(0.25d2+w2/6)3 
along  CB  d^  = d(w-u)/w 

along  EB  d^  = d-u (w-D+d) /w 

approximate  EFB  by  line  EB 


Alternative  case:  for  D-d 

= (w2-0.5d2)/w(0.25d2+w2/6) 
a2  = -12/0. 5d2+w2/3) 
a3  = 4/w(0.5d2+w2/3) 


w;  we  get 


w = D-d 

Take  0 with  respect  to  sector 
AOB  of  radius  (D+d)/4. 

Then  01  = 0 + 4u/(D+d)  rad 


\ 


So  = (0.  ■Su-0.250)[D2-d2  + 2L(D+d)l+L(D-d) 

V = (0. 5n-0. 250)L(D2-d2)0. 5 
o 


Hence 

S = [ 0. 5n-0. 25 (Q+4u/ (D+d) ) !C <D-u) 2- (d+u) 2+2 (L-u)  (D+d ) D+ (L-u)  (D-' 

V =10.25  -0. 125 ( +4u/ (D+d) ) 1 (L-u) [ (D-u) 2- (d+u)2l 

2 3 

and  z = 1-V/V  = a,r  + a_r  + a,r 

o 1 2 3 

for 


ax  = CL(D-d)  + (0.57r-0.25e)  (D+d)  (2L+D-d)  ]/Vq 

a.  =-2f  (^-0/2)  (D+d ) +2L+D-dl/V 
2 o 

a3  = 8/Vo 


(5)  Cylindrical 

(a)  Cord , N=0,  See  2(a) 

(b)  Symmetrical  Perforated  Grains,  Equal  Webs 

(i)  Uninhibited,  Prior  to  Web  Burnthrough,  Any 


s 

o 

V 

o 

Hence  S 

V 

and  z 


0.  25tt  (L-u)  C (D-u) 2-N  (d+u ) 2 1 

1-V/V  = a.r  + a„r2  + a r2 

o 1 2 3 


-2u) 
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Si 


2 (it/  (N-l )-<t>) 


= arcosT 0A/2.0Z  1 
tij  = 2n  / (N- 1 ) - 2$ 

then 

area  of  lateral  surfaces  of  N-l  large,  outer  slivers  S is,  for  any 
• ■»  • k 

regression  distance  w > u > ; 

ST  = (N-l) (L-u) (0.OQ+2X.AQ) 

Li 

area  of  lateral  surfaces  of  N-l  small,  inner  slivers  is,  for  w > u > u^  ; 

= (N-l)  (L-u)  (3  OZ.01) 

area  of  two  end  surfaces  of  N-l  large,  outer  slivers  S£,  for  w > u > u2  ; 

S = 2 (N-l)  (OQ  sin  (0.  50 ) C (2AQ  sin  (0.  5x ) ) 2-  (OQ  sin  (0. 56) ) 2 d*5 

b 

+ 0. 50Q2 (0-sin (6) )-AQ2  (x~sin  (x) ) } 

area  of  two  end  surfaces  of  N-l  small,  inner  slivers  Sg,  for  w > u > u1  ; 

S = 2(N-l)OZ2{sin2(0.5  B.)  1.73206-1.  5CB..-sin  (6, ) 3} 

e iii 

total  area  of  sliver  surfaces  is  therefore 

s = s + s„  + s.  + s 

L E 1 e 


Hence 

V 

= (S„  + S ) (L-u) 0.5 

b e 

z 

= l-v/v 

o 

i.e. 

z 

= 1-0.5(S  +S  ) (L-u)/V 

E 6 O 

s 

A 

c 

where 

V 

o 

= 0.  25tiL  (D2-Nd2 ) as  before, 

similarly  Sq 

Thus 

R 

= S/S 

o 

and 

P 

y 

= dR/dr 
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The  grain  dimensions,  together  with  u = w/2  may  be  used  to  determine  the 
fraction  of  charge  burnt  prior  to  slivering.*  The  grain's  progressivity 
is  optimised  by  increasing  z at  u = w/2,  as  close  to  1.0  as  physically 
possible,  since  sliver  burn  is  a degressive  burn  phenomena.  For  N=7, 
typical  grain  dimensions  are  d=0.1D,  L=2.25D  and  therefore  P =0.37  and 
Z (u=w/2)  = 0.85  with  r=1.532  w/2  at  z=1.0.  * 


Optional  relation  (a)  of  the  general  cubic  type  for  z can  be  used  if  the 
fraction  z^  of  charge  burnt  when  u = w,  and  r^  and  r^  - the  regression 
distances  at  which  the  large  and  small  slivers  are  ail-burnt,  are  known 
i . e . ^ (a) 

2 3 

z = z + 2 =a  + a.r  + a„r  + a,r 
S s o 1 2 3 


where  a^  = 


a 


o 


a = 
o 


0.0 


z — ( 1— z >(w^+2w(r  -w)  J 

w w s 


z - (1-z  ) Cw^  + 2w (r  -w)  1 
w w S 


0 . 5w  • r • r 


a,  = 2 (1-z  ) (r  -2w) 

1 w s 


a = 2 (1-z  ) (r  -2w) 

1 w S 


a.  = 1-z 

2 w 


0.  Sw  < r < r 


1 


r < r < r 
1 2 


0 . 5w  < r < r„ 


Optional 

S 


relation  (b)  During  slivering,  using  a linear  approximation, 


= S Cl-  (z-z  )/(l-z  )3 
w w w 


0.5w  < r < r2 
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generate  a table  of  S versus  z < z < 1.0,  and  interpolate  for  the  required 
value  of  S. 


(i)  r = r when  r = 0.5w  i.e. 

o 

r = 0.125  D-0. 375d 
o 

web  burns  through  to  form  six  large  and  six  small  slivers, 

(ii)  r = r when  8-^  = 0 i.e. 

r = 0. 14433757D-0. 35566243  d 
small  slivers  complete  burning 

(iii)  r = r2  when  8=0  i.e. 

r2  = 0.16930476D-0. 33069524  d 
large  slivers  complete  burning. 
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Similarly,  for  prior  to  sliver  burn 


S = S -z  (S  -S  )/z 
o o w w 


0 < r < 0. 5w 


for  constant  values  S and  z , which  are  calculable  from  the  previous  more 

w w 

accurate  relations.  Option  (b)  enables  significant  saving  of  space  and 
time  for  the  computer  simulation,  and  is  a reasonable  approximation  for 
L/D  >1.0.  It  may  also  be  utilised  for  grains  other  than  7 perforations, 
i . e . any  N . 

(iii)  Uninhibited,  Sliver  Burn,  N=19 


m 


nis 


STAGE  1 


UK 


STAGE  2 


STAGE  3 


STAGE  4 

Slivers  1. 

Inner 

set , 

6 small,  as  for  N=7 

2. 

Second 

. row 

, 6 small,  identical  to 

inner 

set 

3. 

Third 

row, 

12  small,  identical  to 

inner 

set 

4. 

Outer 

set. 

two  stacT'-*  sliver  burn 

(a) 

6 large  slivers 

(b) 

12  residual  slivers 

4 Burning  of  12  residual  slivers,  r^  < r < 


where 


and 


r = 0.0833'D-0.4166‘d 

o 

r 1 = 0.09622504  D - 0.40377496  d 

r2  = 0.10566243  D - 0.39433757  d 

r3  = 0.11286984  D - 0.38713016  d 

Stage  1 - Prior  to  Slivering 

Identical  to  5(b)  (i),  N=19 

Stage  2 

(a)  24  small  slivers 


u < w = (D-5d)/6 


(D-5d)/6  < u < 2r1 


Identical  to  5(b)  (iii),  thus:  S = S.  + S with  (N-l)=24 

s 1 e 

V = 0.5  S (L-u) 
s e 

(b)  6 large  slivers 


OB 

= 

OA  = O.S.  (D-U) 

oc 

= 

CD  = (D+d)/3 

AC 

= 

GE  = 0.5  (d+u) 

CL 

= 

OC/2  = CG 

CH 

= 

OC/4 

OG 

= 

0.866026  CD 

y 

= 

arcosC  (0 C2  + OA2  - 

ac2)/2oc.oa: 

9i 

= 

tt/3-2y 

arc  AB 

= 

OB.  0 x 

OCA 

- 

arcosf (OC2  + AC2  - 

oa2)/2oc.ac: 
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arcos (CH/AC ) 


6 = OCA  - uij  - ti/3 

arc  AE  = AC . 6 


arc  EF  = AC.  (ti  - 2cj^) 

Lateral  area  of  large  slivers  is  therefore  : 


ST  = 6 (L-u)  (AB+2AE+EF) 

JL> 

HE  = AC  sin  (u^) 

MCA  = tt  - OCA 


Area  of  ends  of  large  slivers 


12C0. 5.0A. AB-O. 433013CD2~CG. HE-O. 5GE. EF-AC. AE-CA2MCAl 


SS  - sl  * h1 


vs  =0.5  S£  (L-u) 


Thus  for  stage  2,  with  S and  V identical  to  section  5(b), 

o o 


s = s + s_ 

s S 


V = Vs+VS 

z = 1-V/V 


Stage  3 - Continued  burn  of  6 large  slivers 


where 


= 1-V/V 


i. 


These  are  almost  identical  to  the  large  slivers  of  7 perforated  grains 
as  given  in  section  5(b)  (ii). 


i 

f- 


i 

i 


Hence 


and 


for 


( iv)  Inhibited  External  Wall,  Prior  to  Sliver  Burn 


V = 0.25irL(D2-Nd2) 

o 

S = 0.  5itC  2N  (d+u)  (L-u) +D2-N  (d+u)  21 

V = 0.2571  (L-u)  CD2-N  (d+u)2l 

z = 1-V/V  = a.r  + a_r2  + a r3 

o 1 2 3 

a1  = 2 (2NdL+D2-Nd2)/L(D2-Nd2) 

a2  = 4N(L-2d)/L(D2-Nd2) 

a3  = -8N/L(D2-Nd2) 


(v)  Inhibited  External  Wall  and  Ends,  Prior  to 
Sliver  Burn 

S = irdLN 

o 

V = 0. 25irL(D2-Nd2) 
o 

Hence  S = irL(d+u) 


V 


0. 2 5uL (D2- (d+u) 2N) 


4dN/ (D  -Nd  ) 


Assume  a symmetric  perforation  pattern;  a three-stage  sliver  burn  is 
obtained. 


2 


Stage  1 - Geometry  of  inner  slivers  is  identical  to  case  5(b) (ii) 

for  the  limits  0. 5w  < r < r^. 

S]L  = (N-l)  (L-u)  (3.OZ.0  ) 

and  Se  = 2(N-l)OZ2Csin2(0.5B1)xl.73206-1.5(B1-sin(  ) ) 3 

For  the  huge  outer  sliver  123465,  we  have  : 

Xx  = 2ir/3-w 

a)  = arcos  (AD/AP) 

AD  = 0. 5 OA 

AP  = A Q = 0.  5 (d+u) 

PQ  = AP.Xl 
OA  = 0.25  (D+d) 

Area  of  lateral  surfaces  ST  is  therefore  : 

LX 

SrT  = (L-u)  {ttD+2  (N-l)AP.x,  } 


r 


and  area  of  two  end  surfaces  S is 

El 


EI 


2{  ttD  / 4- 2 (N—  1 ) (OAD)  .-2  (N-l)  (ADP)  -2  (N- 1 ) (APQ) 


sect. 


where  (AOD) , = 0.25  /3  OA^ 


(ADP),  = 0.5  x AD  x AP  x sin  (at) 
A 

2 

(APQ)  = 0.  5 x AP  x v 

sect  A1 


total  area  of  sliver  surfaces  is 


and  since 


where 


S1  + SLI  + Se  + SEI 


Vj  = (Se  + S£I) (L-u) 0. 5 

zl  = 1.0  - (Sg  + S£I) (L-u)0.5/Vq 

V = 0.25irL(D2-Nd2) 
o 


Stage  2 - Continued  burn  of  outer  sliver  under  the  limits  : 
r < r < 0.25 (D+3d) 5w 


where 


and 


S_  = S + S 
2 LI  EI 


V2  = SEI  0,5  (L_U) 


z_  = 1.0  - 0.5  S_T  (L-u) /V 
2 EI  o 


Stage  3 - Residual  burn  of  six  outer  sliver  remnants  under  the  limits: 
w < r < u^  = (5D2  - 6Dd  + Sd2)*5/^ 

3 LI3  LE3 

where  SLI3  = (N-l)  (L-u)  (6 .0Q+2.X- AQ) 

6 = 2 (it/  (N-l)  -<j>) 

<|i  = arcosC  (6A2+0Q2-AQ2)/2.0A.0Q)  3 


OQ  = 0.5  D 
X = OAQ  - at  - ir/3 

OAQ  = arcosC  (OA2+AQ2-OQ2)/2.0A.AQ)3 
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Sink 


* Assume  perforations  are  symmetrically  spaced  and  are  of  equal  diameter. 
For  both  cases  : 

r,  = 0.07735027  d + 0.57735027  w. 

1 l 

D = 3d  + 2w  + 2w . 

o l 

r2  = {0.25(D2-d2)+(l-/3D/2)  (d+wi)2}/(D+(l-/3)d-/3wi) 

Stage  1 


r > 0. 5 w. 

i 


and 


r < (0. 5 w , r ) 
o 1 


cf.  equal  web  case  5 (b) (ii) . 


(a)  six  inner  slivers 


and 


Sx  = (N-l)  (L-u)  (3  OZ.B^ 

S = 2 (N- 1 ) 0Z2{1. 73206  sin2(0.5(3.)  - 1.5  x 
e 1 

x (B1~sin  (6^ ) ) } 


where 


For 


And  if 


Then 

and 
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OZ  = (A Q = AZ  = 0.5) (d+u) 


B1  = 2 ( /(N-l)-^) 
= arcos  (0A/20Z) 
OA  = d + w. 


(b)  large  outer  sliver 
9 = 27I/3-w 


,AD 

to  = arcos  ( /AP) 

AD  = 0.  5 OA 
OA  = d + w. 

l 

AP  = 0.5  (d+u)  = OZ 

S = (L-u)  C tt  (D-u)  + 6. 6.  (d+u)  1 
OJ-» 

PD  = AD.  tan  (to) 

AE  = 2.  AD  = OA 
OD  = OA. 0.5/3 

S^.  = 0.5ir  (D-u)2-3  OA2(/3+tan(io))-(2  -3w)d+u) 
EL 


s = s.  + s + s„  + S„T 

1 e OL  EL 


z = 1.0-0. 5 (S  + S„T)(L-u)/V 
e EL  o 


Stage  2 


(a)  For  r < r < 0. 5 w 
1 o 


S = S„r  + S„ 
OL  EL 


7.  = 1.0  - 0.5  S (L-u)/V 

EL  o 

(b)  If  r.  > 0.5  w and  r.  > r > 0.5  w 
1 o 1 o 

S = S.  + S + ST  + S,, 

1 e L E 

z = 1.0  - 0.5  (S  + S ) (L-u) /V 

e E o 

where  ST  = (N-l)  (L-u) (B.OQ  + 2x-AQ) 

Xi 

AQ  = 0.5  (d+u) 

OQ  = 0.5  (D-u) 

X = OAQ-ui-tt/3 
to  = arcos  (AD/AQ) 

AD  = 0.5  OA 

OA  = d + w. 

l 

OAQ  = arcosC (0A2+AQ2-0G2)/20A.AQD 

B = 2 

<J>  = arcosC  (oa2+oq2-aq2)/2oa.oq: 

and  S£  = 2 (N-l) {OQ. sin (0. 53) C (2AQ. sin (0. 5x) ) 2 - (OQsin  (0. 56) ) 2 Is 

+ 0.5  OQ2  (B-sin  (B) ) - AQ2 (x-sin (x) ) } 


Stage  3 


If  r.  < 0. 5 w and  r„  > r > (0.5  w ,r.) 
1 o2  o l 

5 = SL  + SE 


z = 1.0  - 0.5  S-  (L-u) /V 

r.  O 
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Sliver  Burn,  Outer  Web  Smaller  than  Inner  Web 


r > 0.5 

6 inner  slivers 
6 outer  slivers 


STAGE  2 


j 


(a)  <r 

residual 

slivers 


r 


2 


outer 


r 


2 


SLL  = 6(L"UHVQ  + 2.QR  + TS} 

VQ  = 0. 5 (D-u)B 
B = 2 (ir/6-4.) 


and 


where 

and 


arcosC  (OA2  + OQ2  - AQ2) /20A.0Q1 


OA  = d + w. 

l 

OQ  = 0.5  (D-u) 

AQ  = 0.5  (d+u) 

QR  = 0.  5 (d-u)  (tt-gj) 

(it-oj)  = arcosC (AQ2  + OA2  - OQ2)/2AQ.OAl 

TS  = 7T(d+u)/6 


EE 


0.  5{ttD2-3u  (2D-u)-7tt  (d+u)  2 + 48A*} 


S = STT  + S„„ 

LL  EE 

z = 1.0  - 0.5  S__  (L-u ) /V 
EE  o 


1 2 

A*  = h (d+u)  (d+w^ ) sin  (tt-oj)  + /8(d+u)  uj- 

e = arcosC (0Q2+0A2-AQ2)/20Q.0A3 


Stage  2 

0.5wi  < r < (rirr2) 

c.f.  case  1,  stages  1 and  3(b) 

S = S.  + S + ST  + S„ 

1 e L E 

z = 1. 0-0. 5 (S  +S  ) (L-u)/V 
e E o 


1/8 (D-u)2e 
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